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ABSTRACT 

Double-lined spectroscopic binary systems, containing a Wolf-Rayet and a massive 0-type star, are key objects for the study of 
massive star evolution because these kinds of systems allow the determination of fundamental astrophysical parameters of their 
components. We have performed spectroscopic observations of the star WR68a as part of a dedicated monitoring program of WR 
stars to discover new binary systems. We identified spectral lines of the two components of the system and disentangled the spectra. 
We measured the radial velocities in the separated spectra and determined the orbital solution. We discovered that WR68a is a double- 
lined spectroscopic binary with an orbital period of 5.2207 days, very small or null eccentricity, and inclination ranging between 75 
and 85 deg. We classified the binary components as WN6 and 05.5-6. The WN star is less massive than the O-type star with minimum 
masses of 15±5 Mq and 30±4 Mq, respectively. The equivalent width of the He ii 44686 emission line shows variations with the orbital 
phase, presenting a minimum when the WN star is in front of the system. The light curve constructed from available photometric data 
presents minima in both conjunctions of the system. 
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1. Introduction 

Observational works de dicated to the study of mul t iplicit y 
among massive stars, e.g. lBarba et al.l(l2010ll . lSana et al.l (l2013h . 
and many references there in, provide evidence that the fraction 
of O-hOB systems is high, greater than 50%. This situation is not 
completely replicated by the available data for Wolf-Rayet (WR) 
stars. There are about 50 binary systems detected among the 
galactic WR stars (less than 10% of the known whole sample), 
and only 24 of them are known double-lined spectroscopic sys¬ 
tems (SB2). These systems are key objects because they allow 
the determination of some fundamental astrophysical parame¬ 
ters for their component stars. Moreover, double-lined eclipsing 
systems allow the determination of stellar masses and radii in 
a direct and reliable way, thus providing strong constraints on 
stellar physics and evolutionary models. In this context, the dis¬ 
covery of these kinds of systems is highly relevant. 

In an attempt to search for new binary systems among faint 
WR stars, a spectroscopic monitoring of southern galactic WR 
stars is underway. This survey was started in 2007 and its main 
results published to date are the discovery of the SB2 WR 62a, 
one of the most conspicuous radial-ve locity (RV) variable stars 
in our s ample ([ Collado et al.l 120131 hereafter Paper I), and 
WR 35a (lGamenetal.ll20l4 . 


WR 68a (SMSNP L13, v=14.41) was identihed as a WR star 
bv IShara et al.l (Il999h . No other data or analysis on this star were 
found in the bibliography, thus it was included in our survey. 

In this paper, we present the hrst spectroscopic and photo¬ 
metric analysis of WR 68a, demonstrating that this is an SB2 sys¬ 
tem. The paper is organised as follows: in section 2, we describe 
the observations and data reduction. In section 3, we present the 
spectral analysis of the system’s components, the measurement 
and analysis of radial velocities, and the study of the available 
photometry. In section 4, we summarise our results. 


2. Observations and data reduction 

We obtained thirty-six spectroscopic observations of WR 68a be¬ 
tween 2007 and 2013. The observations were carried out with 
the 2.15-m J. Sahade telescope at the Complejo Astronomico 
El Leoncito (CASLEO, Argentina), the 4-m V. Blanco telescope 
at Cerro Tololo Inter-American Observatory (CTIO, Chile), and 
the 2.5-m du Pont and the 6.5-m Magellan Baade telescopes at 
Las Campanas Observatory (ECO, Chile). See Table[T]for a sum¬ 
mary of the observations. Immediately after or before the stellar 
integration, at the same telescope position, we obtained compari- 
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Table 1. Details of the observing material. 


Date-Obs. 

UT 

n 

Sp. coverage 
[A] 

Observat. 

Telesc. 

[m] 

Spectr." 

Grating 
[1 mm-‘] 

Detector 

Dispersion 
[A pix-‘] 

Resolving ^ 
power 

2007, Mar. 30 

2 

3650-6700 

CTIO 

4 

R-C 

632 

Loral 3k 

1.01 

1300 

2007, Apr. 1-3 

4 

3650-6700 

CTIO 

4 

R-C 

632 

Loral 3k 

1.01 

1300 

2008, Apr. 19-23 

7 

3650-6700 

CTIO 

4 

R-C 

632 

Loral 3k 

1.01 

1300 

2008, Jun. 16 

1 

3650-6750 

LCO 

6.5 

IMACS f/4 

600 

SITe (Mosaic 1) 

0.4 

3300 

2009, Mar. 24-27 

7 

4030-5590 

CASLEO 

2.15 

REOSC SD 

600 

Tekl024 

1.63 

1000 

2009, Jul. 23 

1 

3890-5520 

LCO 

2.5 

B-C 

1200 

Marconi#! 

0.79 

2500 

2009, Jul. 26 

1 

3540-6700 

LCO 

2.5 

B-C 

600 

Marconi#! 

1.55 

1200 

2009, Jul. 27 

1 

3890-5520 

LCO 

2.5 

B-C 

1200 

Marconi#! 

0.79 

2500 

2009, Aug. 14-17 

3 

4040-5700 

CASLEO 

2.15 

REOSC SD 

600 

Tekl024 

1.63 

1000 

2010, Apr. 10-11 

2 

3930-5600 

CASLEO 

2.15 

REOSC SD 

600 

Tekl024 

1.63 

1000 

2010, Aug. 2-8 

4 

3930-5600 

CASLEO 

2.15 

REOSC SD 

600 

Tekl024 

1.63 

1000 

2011, Apr. 6-11 

2 

3930-5600 

CASLEO 

2.15 

REOSC SD 

600 

Tekl024 

1.63 

1000 

2013, May 12 

1 

3930-5600 

CASLEO 

2.15 

REOSC SD 

600 

Tekl024 

1.63 

1000 


Notes, (n): Number of spectra taken per run. Details of the spectrographs can be found in the user manuals of the respective observatories. 
I*"' The spectral resolving power (R=/l/A/l) were measured using Ad as the FWHM of calibration lamp emission lines. 


son lamp spectra. We reduced the spectra using IRAlQ standard introduced by the disentangling method, we could not derive a 
procedures. reliable luminosity class for the secondary spectrum. 


3. Results and discussion 

3.1. The spectrum of WR 68a 

The spectrum of WR 68a shows noticeable emission lines of he¬ 
lium, carbon, and nitrogen ions, and also several absorption lines 
of hydrogen and helium superimposed onto the emissions. As we 
prove below, these absorption lines move in anti-phase with re¬ 
spect to the emissions, indicating that WR68a is a double-lined 
WR -H OB binary system. 

To separate the individual components in the spectrum of 
WR 68a, w e applied a disentangling m ethod similar to that de¬ 
veloped bv iGonzalez & Levatol (12009) . Briefly, the method con¬ 
sists of shifting all the spectra to the common RV for one com¬ 
ponent and co-adding them so that the features of the other 
component are diluted. In a second step, the pure spectrum of 
one component (template) is subtracted to obtain a new RV de¬ 
termination. These steps are repeated in an iterative way. The 
method was only applied to the CTIO and LCO spectra because 
CASLEO data have lower signal-to-noise ratio (S/N), and thus 
they introduce noise into the mean spectra of each component of 
the binary. 

The individual spectra (see Fig. [T]i were used to perform 
a more detailed analysis of both stars. Relative intensities of 
Niii, Niv, Hen, and Civ emission lines in the spectra indicate 
a WN6o_SDectral tvpe (accordin g to the c riteria given bv Smith 
et al.T 996rlT ln good agreement ^thlS hara et al.l 0 9991) . 

Absorption lines of H (Ha, li/3, Hy, H(5), He i (Td4026,4471, 
and 5875), He ii (444542,4686, and 5411), and C iii 45696 emis¬ 
sion were identified in the disentangled spectrum, indicating the 
presence of an O- type star. To determine its spectral type, we 
used the mgb code (iMafz Anellaniz et al.ll2012h . whic h facilitates 
a visual co mparison with standard O-type spectra (ISota et alJ 
l201lLl2014ll . Thus, we classified the O-type star as 05.5-6. Be¬ 
cause of poor S/N in the blue region of the spectrum and prob¬ 
able residuals of the wings of the H5 and He ii 44686 emission 

’ IRAF is distributed by the National Optical Astronomy Observato¬ 
ries, which are operated by the Association of Universities for Research 
in Astronomy, Inc., under cooperative agreement with the National Sci¬ 
ence Foundation. 

^ The “o” means that no hydrogen is observed in the spectrum. 


3.2. Measurement and analysis of radial velocities 


The radial velocities (RVs) were determined by means of cross¬ 
correlation (using the IRAF fxcor task) between the separated 
spectra of each component and the corresponding template. For 
each emission line, an appropriate spectral range was used to 
determine the RV thereof. In the case of the RVs of the absorp¬ 
tion lines, we determined a mean RV using the spectral ranges of 
Hj3, Hy, and H^. The individual RVs of the measured emission 
lines according to the Heliocentric Julian Dates (HJD) of the ob¬ 
servations, including N iv 44058, the blend N iii 444634-40-42, 
Hen 444686 and 5411, and Civ 45812 emissions, and the mean 
RV of H(5, Hy, Hy6 absorption lines are listed in Table |2l 


We searched for periodicities in the RVs of the Hen 44686 
emissionJinej_using_Jhe_jnethodL.dew|^ed_b^Marraco & 
Muzzio (119801) . This algorithm computes, for each trial period, 
the variances from the best-fitting straight line in each phase in¬ 
terval. The most probable period obtained is 5.22 d. Then, by 
means of the gbart0 code, we individually fit orbital solutions 
to the RVs for the emission and absorption lines . As all lines 
show a similar periodicity, we adopted the straight mean of the 
periods obtained from each data set, i.e. P=5.2207+0.0005 d, as 
the binary period. Thus, with the fixed period, we independently 
determined the orbital parameters for the emission and the mean 
of the absorption lines, and we show these in Table[3] The orbital 
solutions for each emission line are illustrated in Fig.|2] 

As expected, each data set presents different systemic ve¬ 
locities, semi-amplitudes, times of periastron passage, and max¬ 
imum RV. For example, the blend N iii 444634-40-42 presents 
the lowest semi-amplitude and the C iv 45812 emission appears 
to have the largest. This is a common behaviour in WN+OB sys¬ 
tems, generally related to the fact that emission lines are formed 
in asymmetric (not barycentr ic) regions at different dep ths in the 
expandingenvelo 2 e_(see_ejgjNiemelu^^o^ 198^Niemela 
et al. ll995h . 


^ GBART is an improved version of the program for the determination 
of the orbital elements for spectroscopic binaries originally written by 
iBertiau & Grobbelil l ll968l) . developed by F. Bareilles, and available at 
http://www.iar.unlp.edu.ar/~fede/pub/gbart. 
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Fig. 1. (a) Continuum rectified spectrum of WR 68a before disentangling; (b) the WR component after removing the 0-type; (c) the remaining 
spectrum of the 0-type star. The spectra are vertically shifted for comparison. 


It is also noteworthy that the RVs of the absorption lines are 
anti-phased with respect to the emission lines, thus indicating 
that they belong to the O-type component of the system. The 
semi-amplitude of its orbit is lower than that of the WN (regard¬ 
less of which emission line is used), which means that the O-type 
star is the more massive star in the system. 

Eccentricities other than zero are not common in binary sys¬ 
tems containing classical He-burning WN component s. Only 
WR 97 is known to have a significant e (0.1+0.04, iGamenI 
(l2004ll l. We checked th e robustness of the foun d e with a test 
of significance given bv iLucv & Sween^(ll971h and it resulted 
to be spurious in all of the emission lines (they do not reach the 
5% level of significance for accepting eccentric orbits), but are 
significant in the absorption lines. This feature should be inter¬ 
preted with tailored models appropriate to this context of close 
binaries, i.e. heating hemispheres by a hot companion, asymmet¬ 
ric emission-line forming regions, colliding wind zone, etc. On 
the other hand, as is shown below, the light curve presents two 
dips at 0.5 phase intervals, also pointing to a circular orbit. In the 
following analysis, a circular orbit is assumed for the WR 68a bi¬ 
nary system. 

To determine the orbital solution of both components simul¬ 
taneously, we assumed that the orbital motion of the WN star 
is better represented by the mean RVs of the nitrogen emission 
lines, and the motion of the O-type star is better depicted by the 
most conspicuous lines in the O-type spectrum. We shifted the 
RVs of emission lines by the difference between the barycen- 
tric velocities obtained for the individual solutions to match the 


Table 4. Cicular orbital elements of WR 68a using the mean of emission 
lines of nitrogen for the WN component and the mean of H/3, Hy, and 
H(5 absorption lines for the O-type component. 


P[d] 

5.2207 (fixed) 

^WR 

295±7 

Kq [km s“*] 

144±9 

Vo Ikm s“*] 

-145±4 

TRVmax Id]* 

5039.56+0.02 

7'wRinfront [d] 

5040.86+0.02 

nwR sin i [Rq] 

33+1 

ao sin i [Rq] 

15+1 

MwRsin^ i [Mq] 

15+5 

Mo sin^ i [M©] 

30+4 

?WR/0 

0.49+0.04 

cr [km s“^] 

32.12 


• HJD-2 455 000 


lines of the the O-type star component. We performed an iter¬ 
ative process to determine the orbital solution, rejecting those 
points whose O-C values were three times larger than the rms of 
the whole solution in each trial. The orbital parameters and the 
RV orbits are shown in Table|4]and Fig.[3j respectively. 

The minimum mass derived for the O-type star is very sim- 
ila r to the theore t ical va lue for an O 5.5-6 V star, as calibrated 
bv iMartins et al.l(l2005h . 31-33 Mq. Therefore, the orbital incli¬ 
nation should be close to 90 deg, and we can expect to observe 
photometric variability (eclipses, for instance). 
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Table 3. Orbital solutions corresponding to the radial velocities of the N iv, N iii, He ii, and C iv emission lines, and Mean Abs. 


Parameter 

Niv 

4057.76 

Niii 

4640.64 

Hen 

4685.68 

Hen 

5411.52 

Civ 

5811.98 

Mean Abs. 

P [days] 



5.2207±0.0005 



Vo [km s^M 

-159±7 

-245±5 

4±3 

22±5 

-373±4 

-146±4 

K [kms-‘] 

314±11 

307±7 

331 ±4 

331 ±7 

352±7 

158±6 

e 

0.05±0.03 

0.08±0.02 

0.05±0.01 

0.04±0.02 

0.06±0.01 

0.20±0.04 

w rdeareesl 

350±37 

206±17 

111±12 

134±30 

157±19 

262±9 

Tperiast [4]* 

5036.9±0.5 

5039.7±0.2 

5038.4±0.2 

5038.8±0.4 

5039.4±0.3 

5038.4±0.1 

TRVmax [d]* 

5037.1 ±0.5 

5036.7±0.2 

5036.8±0.2 

5036.9±0.4 

5037.2±0.3 

5039.4±0.1 

PwRinfront [d]* 

5038.3±0.5 

5038.1±0.2 

5038.1 ±0.2 

5038.2±0.4 

5038.5±0.3 

5035.9±0.1 

a sin i [Rq] 

32±1 

32±1 

34±1 

34±1 

36±1 

16±1 

(T [km s“^] 

31.49 

35.57 

23.51 

37.37 

17.12 

21 


* HJD-2 455 000 



0 0.2 0.4 0.6 0.8 1 1.2 1.4 

Orbital Phase 


Fig. 2. Radial velocities of Niv, Niii, Hen, and Civ emission lines. 
Continuous curves represent the orbital solutions in Table[3] 


3.3. Analysis of the emission-line profiles and available 
photometry 

We analysed the behaviour of the equivalent width (EW) and 
full width at half maximum (FWHM) of some lines at differ¬ 
ent orbital phases. We found that the EW and EWHM of the 
He II /1/14686, 5411 and N iv 44058 emission lines present some 
variability modulated by the orbital period, which is very notice¬ 
able in the He ii 44686 line (see Eig. This line becomes fainter 
around phase 0=0.0, i.e. when the WN star is in the front of the 
system. A similar effect has previously been obser ved in WR 62a 
(see Paper I) and in the intriguing LBV HD 5980 (iFoellmi et al.l 



Fig. 3. WR 68a RV orbit represented by the mean of emission lines of 
nitrogen for the WN component (open circles) and the mean of H/3, Hy, 
and H(5 absorption lines for the secondary component (close circles). 
Points not considered in the solution are plotted with smaller sizes. 



0 0.2 0.4 0.6 0.8 1 1.2 1.4 

Orbital Phase 


Fig. 4. EW and FWHM, as a function of orbital phase, measured in the 
Hen 44686 (close circles), Niv 44058 (open circles), and Heii 45411 
(crosses) emission lines. 


l2008l) . and it has been explained as due to the emission line hav¬ 
ing an additional non-stellar component, which originates in the 
colliding wind region. 

We also analysed th e data collected b y the All Sky Auto¬ 
mated Survey (ASAS) (IPoimahskill200lh . We constructed the 
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Fig. 6. The ASAS V-band photometry of WR 68a. Crosses depict the 
ASAS data, and the circles, the averaged in phase bins of 0.05 with their 
respective standard deviation. 


light curve (see Fig. |6] where we also show the averaged points 
for each bin of <p-Q.Q5), and noted two dips: one in phase 0.0 
(when the WN-type component passes in front of the system) 
and the other, less pronounced, in phase 0.5. To prove it is not 
fortuitous, the ASAS data were searched for periodicities, apply¬ 
ing the same algorithm used for the RVs. The most likely period 
found is P- 5.22 d, which is similar to the spectroscopic period 
(both periodograms are compared in Fig.|5]). 

The similarity between the light curves of WR62a and 
WR 68a is notorious, however, the dip at 0=0.0 in the light curve 
of WR 68a is shallower and wider than that observed in WR 62a. 
The resemblance of these light curves to that of V444 Cyg, a 
WN5H-06 binarv svstem, is also remarkable (see Figure 2 of Eris 
& Ekmekci 1201 ID. The photometric data of the ASAS survey 
are not suitable for a detailed analysis because of the faintness 
of WR 68a. Notwithsta nding that, we adjusted a si mple Wilson- 
Devinney (WD) model (IW ilson & Devinnevlll971h by means of 
the PHOEBE code (iPrsa & Zwitted 2005 ). to the RV of both stars 
and the photometric data together. We adopted t he theoretical 
stellar mass and radii for an O 5.5-6 V star from iMartins et al] 
(l2005h and adjusted the stellar parameters of the WN6 compo¬ 
nent. We found that the light curve can be reasonably htted with 
an orbital inclination ranging between 75 and 85 deg. These val¬ 
ues imply for the WN component a radius between 4.5 and 7.4 
Rq, compa rable with those e s timate d for the WN5 component of 
V444 Cyg (lEris & Ekmekcil (1201 Ih and references therein). 


WR 68a becomes in an important astronomical target, as 
accurate multi-band photometric and higher resolution spectro¬ 
scopic follow-up observations will allow to determine the stellar 
parameters, as radius and masses, for both massive components. 
WR 68a is thus a target of choice to obtain direct properties of 
stars in a WR-hO system, hence, to obtain valuable constraints to 
confront evolutionary models of massive stars. 
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4. Conclusions 

We have discovered that WR68a is a double-lined binary sys¬ 
tem composed of a WN6 star and an O 5.5-6 type component 
in a orbital period of 5.2207 d. We found that the O-type star is 
more massive than the WN component with minimum masses of 
30 + 4 Mq and 15 + 5 Mq for the O-type and WN stars, respec¬ 
tively. 

We also detected a minimum in the EW of the He ii 44686 
emission when the WN star is in front of the system, likely due 
to the eclipse of an additional component of this line originated 
in the colliding winds region. From the analysis of the available 
photometric data (ASAS), we found a variability, i.e. two min¬ 
ima in the light curve during both conjunctions of the system, 
which we interpreted as eclipses. 
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Fig. 5. Periodograms obtained for the ASAS data (top) and for the 
Hen 44686 emission line (bottom). 
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Table 2. Observed heliocentric radial velocities of WR 68a. 


HJD 

Phase 

Niv 

Niii 

He II 

Hen 

Civ 

Mean Abs. 

2450000H- 


4057.76A 

4640.64A 

4685.68A 

5411.52A 

5811.98A 

H(5-Hr-H/3 



[km s"'] 

[kms-'] 

[kms-'] 

[km s"'] 

[kms-'] 

[kms-'] 

4189.83 

0.99 

-238+2 

-173+2 

58+1 

46+5 

-471+3 

-113+3 

4189.85 

0.99 

-239+1 

-203+2 

41+2 

30+1 

-491+4 

-158+3 

4191.86 

0.38 

54+2 

-85+2 

190+1 

274+2 

-68+5 

-290+6 

4192.89 

0.58 


-450+1 

-268+4 

-162+7 

-441+3 

-110+6 

4193.67 

0.72 

-425+1 

-550+2 

-309+3 

-301+3 

-702+3 

10+9 

4193.90 

0.77 

-483+2 

-542+1 

-297+1 

-327+4 

-742+3 

-15+9 

4575.90 

0.94 

-297+1 

-281+2 

-64+1 

-113+3 

-571+3 

-36+3 

4576.86 

0.12 

45+2 

-17+1 

256+1 

243+1 

-213+3 

-209+13 

4577.75 

0.29 

150+2 

-38+2 

272+1 

340+1 

-32+5 

-283+17 

4578.66 

0.47 


-299+3 

6+1 

103+1 

-230+2 

-268+7 

4578.89 

0.51 

-188+4 

-365+3 

-157+2 

-70+6 

-249+2 

-161+5 

4579.61 

0.65 

-365+4 

-545+2 

-303+1 

-239+1 

-587+2 

-32+10 

4579.88 

0.70 

-478+5 

-592+1 

-342+1 

-306+2 

-693+3 

-2+14 

4633.69 

0.01 

-213+3 

-160+3 

66+3 

121+3 

-411+7 

-113+21 

4914.87 

0.87 


-415+3 

-177+3 

-270+2 



4915.85 

0.06 


18+1 

202+1 

231+3 



4915.88 

0.06 


-28+1 

214+1 

194+4 



4916.81 

0.24 


-8+3 

314+1 

308+1 



4916.84 

0.24 


33+4 

329+2 

326+1 



4917.80 

0.43 

66+2 

-114+3 

120+1 

157+1 


-304+13 

4917.83 

0.43 

49+2 

-123+1 

111+2 

148+1 


-283+12 

5035.65 

0.00 

-256+2 

-128+2 

93+1 

66+2 


-103+5 

5038.62 

0.57 

-228+5 

-392+1 

-230+4 

-225+5 

-399+2 

-76+9 

5039.64 

0.77 

-444+1 

-531 + 1 

-283+1 

-338+3 


28+20 

5059.56 

0.58 


-440+3 

-194+3 

-131+3 



5060.53 

0.77 



-304+2 

-296+2 



5060.56 

0.77 



-317+1 

-234+1 



5296.85 

0.03 


-90+2 

135+2 




5297.81 

0.22 



317+1 

367+1 



5410.52 

0.81 


-557+2 

-239+1 




5411.55 

0.00 


-158+2 

128+4 




5413.50 

0.38 



197+1 

258+1 



5416.51 

0.95 


-156+2 

-30+2 




5657.88 

0.19 


73+2 

309+3 

268+4 



5662.74 

0.12 



284+2 

216+1 



6425.80 

0.28 


2+3 

322+1 

353+3 



FWHM^ [A] 


14.6 

30.0 

22.5 

21.6 

34.5 


EW^ [A] 


-4.6 

-8.5 

-31.4 

-5.3 

-5.6 



Notes: Radial velocities errors listed in Table are provided by the fxcor task of IRAF. The errors are 
calculated from the asymmetric noise of the cross-correlation function and of the fitted height of the 
peak. 

("i The mean FWHM and EW of each emission line were measured in the disentangled spectra. 











